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ABSTRACT 

We present broad band analysis of the Be/X-ray transient 1 A 1 118-61 by Suzaku at the peak of its 3rd ob- 
served outburst in January 2009 and 2 weeks later when the source flux had decayed by an order of magnitude. 

The continuum was modeled with a cutoff pi model as well as a compTT model, with both cases requiring 
an additional black body component at lower energies. We confirm the detection of a cyclotron line at ~ 55 keV 
and discuss the possibility of a first harmonic at ~ 1 lOkeV. Pulse profile comparisons show a change in the 
profile structure at lower energies, an indication for possible changes in the accretion geometry. Phase resolved 
spectroscopy in the outburst data show a change in the continuum throughout the pulse period. The decrease in 
the CRSF centroid energy also indicates that the viewing angle on the accretion column is changing throughout 
the pulse period. 

Subject headings :X- rays: stars — X-rays: binaries — stars: pulsars: individual (1A 1 1 18-61) — stars: mag- 
netic fields 


1. INTRODUCTION 

The Be/X-ray binary transient 1A 1 1 18-61 was serendipi- 
tously discovered during an observation of the nearby binary 
system CenX-3, when an outburst was detected in Decem- 
ber of 1974 by the Ariel-5 satellite (Eyles et al. 1975). A 
second, similar outburst occured in January of 1992 and was 
observed by the Burst and Transient Source Experiment on 
the Compton Gamma Ray observatory CGRO/BATSE (Coe 
et al. 1994). The measured peak flux was ~ 150 mCrab for 
the 20-100keV energy range, similar to the 1974 outburst. 
The source showed an elevated emission throughout the next 
~ 30 days after the decay of the main outburst (see Coe et al. 
1994, Fig. 1). The third and most recent outburst occurred on 
2009, January 4 and was detected by the Swift Burst Alert 
Telescope BAT (Mangano et al. 2009; Mangano 2009). It 
reached a peak flux of ~ 500 mCrab in the 1 5-50 keV energy 
band. This last outburst was monitored with Swift and the 
Rossi X-ray Timing Explorer ( RXTE ) as well as with two 
long Suzaku pointings and one observation with INTEGRAL 
during a flaring episode ~ 50 days after the peak of the main 
outburst (Leyder, Walter & Lubinski 2009). 
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Pulsations with a period of 405.3 ±0.6 s were observed dur- 
ing the 1974 outburst and were initially attributed to the or- 
bital period of two compact objects (Ives, Sanford & Bell 
Burnell 1975). Shortly afterwards it was suggested that the 
period stems from a slow rotation of the neutron star (NS) it- 
self (Fabian 1975). During the 1992 outburst pulsations with 
a period of ~ 406.5 s were detected up to 100 keV, show- 
ing a broad, asymmetric, single peak pulse profile above 
the lowest BATSE energy of 20keV. The pulse period de- 
creased throughout the decline of the outburst with a rate of 
~ -0.016s/day and it appeared constant at ~ 406.35 s for the 
time of the elevated emission. During the 2009 outburst a sim- 
ilar period evolution was observed with RXTE resulting in a 
pulse period of P sp i n = 407.719s, and P sp ; n = -4.6 x 10 -7 s/s 
« -0.04s/day (Doroshenko et al. 2010). Furthermore, the 
lower energies showed a more complex pulse profile with two 
peaks below an energy of ~ 10 keV. Due to the short duration 
of the Suzaku observations with respect to the pulse period, 
the derived RXTE values are used for determining the pulse 
profile and phase resolved spectra in this paper. 

The optical counterpart was identified as the Be-star Hen 
3-640/Wray 793 by Chevalier & Ilovaisky (1975) and clas- 
sified as an 09.5IV-Ve star with strong Balmer emission 
lines indicating an extended envelope by Janot-Paeheco, Ilo- 
vaisky & Chevalier (1981). The overall spectrum was found 
to be similar to other known Be/X-ray transients, such as 
X-Per and A 0535+26 (Villada, Giovannelli & Polcaro 1992, 
and references therein). The distance was estimated to be 
5±2kpc (Janot-Pacheco, Ilovaisky & Chevalier 1981) and 
was confirmed by Coe & Payne (1985), along with the spec- 
tral type classification, using UV observations of the source. 
EXOSAT observed X-ray emission from 1A 11 18-61 be- 
tween outbursts (Motch et al. 1988), thus indicating a con- 
tinuous low level of accretion. Rutledge et al. (2007) reported 
on pulsations in the low luminosity state observed with Chan- 
dra, making it only the third known HMXB transient after 
A 0535+26 and 4U 1145-619 for which this behavior has been 
observed. 

A study of the H a emission line before and during the 1994 
outburst (Coe et al. 1994) showed a strong correlation be- 
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FIG. 1.— Swift /BAT lightcurve during the 1A 1 118-61 outburst. The 
dashed lines mark the Suzaku observation as well as 2009 Jan 1 . 

tween its strength and the observed X-ray activity, indicating 
the existence of a very large disk around the star. The fluctu- 
ation of the equivalent width of H a indicated possible insta- 
bilities in the disk, which were enhanced when the star passes 
through the periastrion. The analysis of the UV continuum 
and line spectra has indicated that the photospheric emission 
from the Be star was not affected by the X-ray radiation, sim- 
ilar to the case of A 0535+26 (de Loore et al. 1984). 

Until recently, the orbital period of 1A 1118-61 was not 
measured and assumed values were of the order of 350 days, 
based on the P S pi„-P 0 rb relation (Corbet 1986), and ~ 
585 days based on the equivalent width of the H a line (Reig, 
Fabregat & Coe 1997). Staubert et al. (2010) analyzed the 
pulse arrival time in RXTE data throughout the 2009 outburst 
and established an orbital period of ~ 24 days with a very cir- 
cular orbit around the Be-companion. 

Spectral fitting during the 1974 outburst indicated a variable 
power law index, where the lowest value was ~ 0.9 during the 
peak of the outburst, and ~ 1.1 before and after (Ives, Sanford 
& Bell Burnell 1975). Coe et al. (1994) described the 1992 
CGRO/BATSE data with a single-temperature, optically- 
thin, thermal bremsstrahlung (OTTB) with a temperature of 
15 keV. Due to the different energy ranges of the instruments 
and the fact that pulse-off-pulse data have been used, these 
results can not be directly compared with the 1974 observa- 
tion. 

In the 2009 outburst Doroshenko et al. (2010) discovered 
a cyclotron resonance scattering feature (CRSF) at ~ 55 keV. 
CRSFs can be used to deduce the strength of the magnetic 
field at the pole region of the NS and have been observed in 
multiple sources at energies from ~ 1 5 keV up to ~ 50 keV 
(Cobum et al. 2002). 

Suzaku observed 1A 1 118-61 twice in 2009, once during 
the peak of the outburst and again ~ 13 days later when the 
flux had returned to a level slightly above the quiescent state. 
In this paper we will present the first detailed 0.7-200 keV 
broad band spectrum for this source, confirming the detec- 
tion of a CRSF at ~ 55 keV. The paper is divided as fol- 
lows: Section 2 describes the data accumulation and extrac- 
tion; Section 3 presents the broad band phase averaged spectra 
for the outburst (0.7-200 keV) and for the second observation 
(0.7-70 keV), and discusses the CRSF and Fe-Line complex; 
Section 4 compares energy resolved pulse profiles for both 
observations and discusses the phase resolved analysis for the 
outburst data; finally Sections 5 and 6 present a discussion of 
our observations and summarize our findings. 


2. OBSERVATION AND DATA REDUCTION 

A sudden increase in activity of 1A 1118-61 was de- 
tected with the Swift/BAT instrument on January 4th, 2009 
(Mangano et al. 2009). The count rate increased steadily 
until January 15th when it reached the maximum value of 
~ 500mCrab and then decayed until January 27th where it 
showed a low emission level with periods of flaring (Leyder, 
Walter & Lubinski 2009) until mid March and then returned to 
quiescence (see Fig. 1). Suzaku observed 1 A 1 118-61 during 
the peak of the outburst on January 15th, 2009 (MJD 54846.5, 
ObsID 403049010) with both of its main instruments: the X- 
ray Imaging Spectrometer (XIS; Mitsuda et al. 2007) and the 
Hard X-ray Detector (HXD; Takahashi et al. 2007). A sec- 
ond observation was performed on January 28th, 2009 (MJD 
54859.2, ObsID 403050010), ~ 13 days after the main out- 
burst, at the beginning of the period of elevated emission. 

Both observations were performed using the HXD nomi- 
nal pointing to minimize the pile-up fraction in the XIS in- 
struments and to enhance the HXD sensitivity for a possible 
CRSF detection. Of the four original XIS instruments only 
XIS 0, 1 and 3 were functional during the observing time. 
XIS 0 and 3 are both front illuminated (FI) CCDs and XIS 1 
is a back illuminated (BI) CCD. To reduce pile-up, the XIS 
instruments were operated with the 1/4 window option with a 
readout time of 2 s and the burst option with only 1 s accumu- 
lation time for each readout cycle, reducing the exposure time 
for the XIS instruments by 50%. The data were taken in the 
3x3 and 5x5 editing modes which were then combined for 
the final spectral analysis. 

For the extraction the Suzaku FTOOLS version 16 (part 
of HEASOFT 6.9) was used. The unfiltered XIS data were 
reprocessed with caldb20090402 and then screened with the 
standard selection criteria as described in the ABC guide 11 . 
Each detector and editing mode combination was extracted 
independently and individual response matrices and effective 
area files were created. For the final spectra the data from 
both FI detectors were combined (XIS 0 and 3) and the re- 
sponse matrices and effective areas were weighted according 
to the accumulated exposure time of the different modes. The 
XIS data were grouped so that the minimum number of chan- 
nels per energy bin corresponded to at least the half width half 
maximum of the spectral resolution, i.e. grouped by 8, 12, 14, 
16, 18, 20, 22 channels starting at 0.5, 1, 2, 3, 4, 5, 6, and 
7 keV, respectively (Nowak, priv. com). 

For the HXD data the Suzaku team provided the tuned PIN 
non X-ray background 12 (NXB). Following the ABC Guide 
the cosmic X-ray background (CXB) was simulated and the 
exposure time was adjusted in both backgrounds by the pre- 
scribed factor of 10. The PIN data were grouped so that at 
least 100 events were detected in each spectral bin. 

GSO data were extracted using the FTOOL hxdpi with 
the newest gain calibration file from April, 2010. The NXB 
background files version 2.4 created by the Suzaku HXD in- 
strument team were used. The data were then binned to 64 
bins according to the Suzaku homepage 13 . GSO data in the 
range 70-200 keV were used in the spectral analysis. 

These selection criteria resulted in ~ 25 ks exposure time 
for each XIS instrument and ~ 49 ks for the HXD instruments 
in the first observation. The second observation had an expo- 
sure of ~ 21 ks and ~ 29 ks for the XIS and HXD instruments, 

! 1 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

u http://heasarc.nasa.gov/docs/suzaku/analysxs/pmbgd.html 

1 J http://heasarc ,gsfc .nasa.gov/docs/ suzaku/ analysis/gsojnewgain .html 
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respectively. 

2.1. Pileup correction 

For bright sources, such as 1A 11 18-61 
during the outburst, a strong pileup is ex- 
pected. Following the description provided at 

http: // space. mit . edu/ASC/ software/suzaku/ 
the S-lang routine aeattcor . si was used to improve the 
attitude correction file and the point spread function of the 
events. Then the tool pile_estimate . si was applied 
to produce a 2 dimensional map of the pileup fraction. The 
maximum values for pileup fractions were 10% and 15% 
for XISO, 15% and 16% for XIS1, and 18% and 21% for 
XIS3 for the 3 x 3 and 5x5 editing modes, respectively. 
Regions with a pileup fraction above 5% were excluded 
from the extraction for each individual source event file. 
For the second observation the calculated maximum pileup 
fractions were < 5% and no regions had to be excluded for 
the extraction. 

3. PHASE AVERAGED SPECTRUM 

For the outburst observation we extracted broad band XIS 
spectra for 0.7-12keV (0.7-10keV for the BI XIS1), PIN 
spectra for 12-70keV and GSO spectra for 70-200 keV. For 
the second observation the GSO spectrum was not well con- 
strained and therefore was not included in the analysis. The 
final model included the Galactic and intrinsic absorption, a 
continuum that steepened at higher energies, an iron line com- 
plex, and a CRSF. Excessive low energy flux was modeled 
with a black body component. In addition a lOkeV absorp- 
tion feature was required in the outburst observation. 

The Galactic column density was modeled with a single 
photon absorption (phabs) component, where the column 
density was confined between 1.1 and 1.4 x 10 22 cm -2 . The 
value determined by the NASA Ah Tool 14 for 1A 11 18-61 
is 1.22 x 10 22 cm~ 2 . The IVh value was left free in the fits 
for the first observation. For the second observation the val- 
ues showed larger error bars with the cutoff pi model and 
had to be fixed when using a compTT model. The intrin- 
sic column density was modeled with the partially covered 
photon absorption model pcf abs to take the flux at lower 
energies (<1 keV) into account. All modeling was performed 
using the wilm abundances (Wilms, Allen & McCray 2000) 
with the vern cross-sections (Vemer et al. 1996). The con- 
tinuum was modeled using a power law with an exponen- 
tial cutoff (cutof fpl). Using a power law with a Fermi- 
Dirac (FD) cutoff, one of the empirical models often applied 
to accretion powered pulsars, resulted in a cutoff energy of 
~ 0, making the FD cutoff effectively a cutof fpl, where 
the variable E cut actually reflects the folding energy of the 
model. The additional lOkeV feature was modeled using a 
negative, broad (a = 1 -2keV) Gaussian component. Such 
a feature has been previously observed in different sources 
and it is believed to stem from an improper modeling of the 
continuum (see §6.4 in Cobum et al. 2002). More recent ex- 
amples of such a feature are, e.g. 4U1907+09 (Rivers et al. 
2010) and Cen X-3 (Suchy et al. 2008). Although this fea- 
ture sometimes is described as a broad emission line, e.g. 
Suchy et al. (2008), in this case a negative Gaussian line at 
~ 10 keV fits best. When including this line in the overall best 
fit, the y 2 /dofs (degrees of freedom) decreased from 1 173 / 
477 to the best fit value of 678 / 474, respectively. Due to 

14 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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Fig. 2.— Broad band phase averaged spectrum of the outburst observation 
using the cutof fpl model. The residuals are shown without CRSF (b), 
with one CRSF at 55keV (c) and with the second CRSF at 112keV. The 
Inset shows the Fe line region and the best fit model for XIS1 and XIS 03. 

the small differences in the instrument response, all detectors 
were coupled with a normalization constant, which was set to 
1 for the combined XIS 03 detectors. The final model is of 
the form const*phabs*pcf abs* ( cutof fpl*gabs+ 
2*Gauss Fe )+Gaussiokev+3*Gausscai. The best fit pa- 
rameters are mentioned in the text and partially summarized 
in Table 1 . 

A second approach for the continuum was initiated by 
Doroshenko et al. (2010) using the Comptonization model 
compTT developed by Titarchuk & Hua (1995) instead of 
the cutof fpl model. The parameters of the compTT were 
again left independent between the FI XIS03 and the BI XIS 1 , 
but only small differences could be observed. The negative 
Gaussian at lOkeV was not necessary for this model. When 
comparing the spectral parameters in the outburst between 
these two models, i.e. CRSF, Ah and the spectral lines, to the 
previously described cutof fpl, the values were consistent. 

A comparison of the compTT values, such as the photon 
seed temperature To , the plasma Temperature kT c or the op- 
tical depth t p showed consistency between the spectra of the 
outburst as seen by Suzaku and RXTE (Doroshenko et al. 
2010 ). 

When comparing the outburst spectrum with the second ob- 
servation, significant changes in the continuum parameters are 
observed. In the cutof fpl model the power law index soft- 
ened from 0.350^;{g 2 to 1 .0 The cutoff energy £ cut in- 
creased from 1 8.7ig;| keV during the outburst to 24.7^lJ;| keV 
in the second observation. The temperature of the additional 
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TABLE 1 

Phase averaged spectral parameter. Model described in the text. 


Parameter 

phabs .VHflO^/cm 2 ] 
pcfabs WH[10 22 /cm 2 ] 
covering fract. 
blackbody kT [keV] 
blackbody norm [10 -3 ] 
£cut [keV] 

r 

A p) [lO- 2 ]t 


Outburst” 
1 74 +o.oi 

1 1 9+^.2 

066^> 
u - oo ~0.01 
0 «1+0;0i 

18.7; 
0.350; 
9.55 


Mi 

ksr 


£crsf [keV] 

(tcrsf [keV] 
tcrsf 

E Pe K a [keV] 

EfcKg [keV] 

Eq. Width K Q / Kg [eV] 
Fluxa-iokevtlO 36 ergs/sec ]* 
Cxisi/Cpin/Coso 
X 2 /dofs 


58.2^ 

W.Hj 

60.1«| 

A 414+0 006 
7 -t o+U.03 

5l!t/8! 3 

8.79 

0.98/1.10/0.82 

678/474 


Cutoffpl 
2nd Obs. 
1.34 frozen 
13 5 40 - 4 

Q -71+0.01 

0 1 5+^-Q 1 

u - 1 J ~p.oi 

7 7+0.2 

1 02 +<M 

fo$$ 4 


58.2 frozen 
14.1 frozen 
15.4^ 

6 49+0.O1 

7.13 frozen 

46!j§/ 


Parameter 


0 /0+2 
0/ -2 
1.72 

1 . 02 / 1 . 01 /- 

574/431 


2nd Obs, free CRSF 

1-43111 


CompTT 


13.6 

Q 71+ffOl 

0^:81 

24.7+H 

1 or™-® 

skf 


4 7 - 4 fl 

6 0 + ^ 
64 2 *.ii 

O, ^ Z -0.0i 

7.13 frozen 
46!}«/8!l 
1.72 

1 . 02 / 1 . 01 /- 

567/428 


compTT To [keV] 
compTT kT [keV] 
compTT t p 
compTT norm ! 


Outburst 
1 7Q+0.03 
1 ■'^“-0.02 
12.00™- 4 
0 71+® 

0 181 
U,JJ -0.01 
12.85i 

1 34 + 0-04 

7 > 5 

6 00*^ 
°- uu ~0.18 

0.096^;1 3 

54.5jT 

10-34? 

23.6 1> 0 „ 4 
6 40+0:01 

,75*84 

/ ' i +-0.03 

51^/10t 3 

9.15 

0.99/1.18/1.10 
680 / 474 


2nd Obs. 

1 .28 frozen 
15 l"* 0 - 1 
0.76^»> 
0.17" 
2.89" 
n 44^ 

u.uuz^QQi 

54.5 frozen 
10.3 frozen 
77 i+3.2 

6 42^1 

7.13 frozen 
49!jg/ll+« 
1.72 

0.95/ 1.00/- 
567/428 


NOTE. — a Best fit values including first harmonic at ~ 1 lOkeV, see text for details, f Units in Photons keV 1 cm 2 s _1 , f unabsorbed flux using a distance of 5 kpc 


-r-T 
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FIG . 3 . — Broad band phase averaged spectrum of the second observation 
using the cutoffpl model. Only XIS and PIN data were used. Residuals 
without a CRSF (b) and with one CRSF at ~ 47 keV. The Inset shows the Fe 
line region and the best fit model for XIS1 and XIS 03. 

black body decreased from 0.51^;^ 3 keV to 0.1 5i§;§} keV 
after the outburst. The intrinsic column density and cov- 
ering fraction increased both from H.2^' 2 x 10 22 ctrf 2 and 
0.66^oi in the outburst to 13.6ig;| x 10 22 cnr 2 and 0.71^ 
in the second observation. 

In the compTT model the changes between the two obser- 
vations were similar. The intrinsic column density and black 
body values were consistent with the values obtained with 


the cutoffpl model, where the intrinsic column density is 
12.(Sgj x 10 22 cm' 2 and 15.1^ x 10 22 cnT 2 and the cover- 
ing fraction is 0.7 l^J and 0.76^;oJ for the outburst and the 
second observation, respectively. When establishing the er- 
rors for the other spectral parameters of the continuum, these 
values had to be frozen to avoid a drift into an unreasonable 
part of the parameter space. The black body temperature de- 
clined from 0.53 ^-qI keV to 0. 171 qq 5 keV, values as in 
the cutoffpl model. The compTT values changed signif- 
icantly between the two observations. The photon seed tem- 
perature decreased from 1 ,34to'o5 keV t0 0-4+qmkeV. The 
electron plasma temperature tcT e increased from 7.61{j3keV 
to 182^9g 8 keV, whereas the optical depth of the plasma r p de- 
creased from 6.00^;}| to 0. 68^26. This behavior reflects the 
known negative kT e and t p correlation (see Section 5.2). 

3.1. Fe line component 

An Fe K Q emission component was detected at 
6.414^Qg|keV in the outburst data (see Fig. 2, Inset), 
as well as at 6.421^} keV in the 2nd observation (see Fig. 3, 
Inset). For the best fit values the centroid energies of XIS03 
and XIS 1 were left independent, due to a known shift in the 
energy calibration. Although the centroid energies for XIS1 
generally show a slightly lower (30 eV) energy than the other 
XIS instruments 15 , our XIS 1 data instrument actually showed 
a 20-30 eV higher centroid energy. Residuals at ~ 7 keV 
were observed and were modeled with an additional Gaussian 
line with a best-fit energy of 7.13 ± 0.03 keV, indicating the 
existence of an Fe K g line. The width a of the Fe Kg line was 
set to the Fe K a line value, and the normalization was left in- 
dependent. The obtained value for the Kg normalization was 
consistent with the expected 1 2% of the K Q normalization. In 
the second observation the second line was very weak and its 
parameters could not be properly constrained. In this case the 
line centroid was frozen at the outburst value of 7 . 1 3 ke V. The 
width was again coupled with the Fe K Q value. The residuals 
decreased from \ 2 = 721 for 477 dofs to the overall best 

15 http://www.astro.isas.jaxa.jp/suzaku/process/caveats/caveats_xrtxis06.html 
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fit values of 678 for 474 dofs, including two CRSFs, when 
the Fe Kg line was included and an F-test showed only a 
2.8 x 10 -6 probability that this improvement is of a statistical 
nature. 16 Another possibility for modeling the observed 
residuals at ~ 7 keV is the addition of an Fe K-edge, located 
at 7.1 keV. When including the edge, the improvement is only 
marginal (x 2 of 712 for 474 dofs) and the best fit energy is 
~ 7.6 keV. Freezing the energy to the value of 7.1 keV did 
not improve the fit. The remaining residuals indicated that an 
additional Fe Kg line would still be necessary. The measured 
energies of the Fe K Q and Kg lines were slightly higher than 
the expected laboratory values, but still in the same order of 
magnitude of the energy scale uncertainties of the instrument 
( ~ ± 20 eV @ 6 keV according to the ABC guide). 17 

Comparing the equivalent widths (EQWs) of the Fe K Q 
and Kg lines in the two observations showed no significant 
change during the outburst. The Fe K Q line showed an EQW 
of 5 1 ± 4 eV for both models, whereas the Fe Kg line showed 
an EQW of 8 ±3 eV for the cutoff pi model and 10±3eV 
in the compTT model. For the second observation, the Fe K Q 
EQWs were 46 ± 10 eV and 49 ± 10 eV, and the Kg EQWs 
were 8±2eV and 11 ±9 eV for the cutoff pi and compTT 
models, respectively. The similarity of the Fe line EQW be- 
tween both observations indicates that the source of the Fe 
lines is near to the source of the continuum (see discussion). 

3.2. Cyclotron resonance scattering feature 

In the burst observation a strong residual at ~ 50 keV was 
observed in the HXD (Fig. 2, b). When including a CRSF-like 
feature with a Gaussian optical depth profile (gabs) the best 
fit was obtained with a centroid energy of 58.27$;$ keV, signif- 
icantly improving the x 2 from 1728 with 480 dofs to 752 for 
477 dofs for the cutoff pi model, confirming the discovery 
of the CRSF in the RXTE data by Doroshenko et al. (2010). 
The width of the line was 14. llH|;f keV and the optical depth 
r was 60. llf |. Using the compTT model, similar values for 
E c yc and a were obtained for the outburst data. In this case the 
optical depth was found to be 23.67g°9 4 . 

In the cutoff pi model, a second absorption line, using 
the gabs model, in the GSO data of the outburst spectrum at 
~ 1 10 keV improved the x 2 to 678 with 474 dofs (see Fig. 2c). 
An F-test probability of 0.144 indicates that the line may not 
be significantly detected. However, an increase of the GSO 
Non X-ray Background (NXB) by 1 or 2% decreased the GSO 
counts to a non-detectable value. Another possible system- 
atic effect is the decay of the 153 Gd instrumental line, which 
is activated when the satellite passes through the south at- 
lantic anomaly (SAA), and thus creating a background line at 
~ 1 50 keV. This could impose a deviation at ~ 110 ke V if not 
accurately represented in the modeled background. Therefore 
caution is advised when interpreting this feature. 

In the lower luminosity observation only an indication of 
the fundamental CRSF was observed (see Fig. 3). The line 
was fitted in 2 steps, where first the centroid energy and width 
were fixed to the outburst values, and the optical depth was 
left as a free parameter. With the addition of this CRSF, the 
cutoff pi model improved slightly from a x 2 /dofs of 605 

16 See however Protassov et al. (2002) about the usage of the F-test in 
line-like features. 

17 An inflight energy calibration using the Mn-calibration sources of the 
XIS detectors could not be performed, due to the usage of the 1/4 window 
mode, which excludes the comer regions where the calibration sources are 
located. 


/ 432 to 574 / 431 with an optical depth of r = 15.47$;$$. In 
the second step, the CRSF parameters were left free and the 
centroid energy decreased to 47.47$;$ keV, whereas the opti- 
cal depth decreased to r = 6 . 07 $$ and the width decreased 
to cr = 5.7*j ®keV. With a x 2 /dofs of 567 / 428 the best fit 
improved further (see Tab. 1 for details). To test the signif- 
icance of the CRSF in the second observation, 100000 PIN 
spectra were Monte Carlo simulated using the null hypothesis 
approach where a spectrum, including Poisson noise, was cre- 
ated using the best fit parameters without an CRSF line. The 
simulated spectra were then fitted with a continuum model 
and an additional gabs absorption line to test how often such 
a line could be detected in the spectral noise. The width was 
constrained to a value of a between 3.5 keV and 8.5 keV, so 
that neither very small features nor broad parts of the con- 
tinuum were modeled. Out of the 100000 simulations, 42% 
of the best fit results showed a non-zero value of r. When 
comparing the simulated r distribution with the best fit result 
of r = 6.07$$, ~ 96% of all simulated fits showed a value be- 
tween 0 and the measured value. When taking the errors of the 
measured data into account (r < 3.9) the number is reduced 
to~90%. 

For the compTT model, the addition of the fundamental 
CRSF with the frozen energy and width created residuals in 
the 1 0 - 20 keV range , which could be significantly reduced by 
decoupling the optical depth of the compTT model in the PIN 
data from the XIS data. The best fit optical depth of the CRSF 
line was 27. changing the x 2 / dof values from 639 / 429 
to 567 / 428, respectively. When leaving the CRSF parame- 
ters free, the width of the CRSF increased to > 30keV, mak- 
ing it part of the continuum and tampering with the results. 
Freezing the CRSF parameters to the values obtained with the 
cutoff pi model did not result in a satisfactory overall fit. 

An alternative approach to describe the CRSF is to replace 
the gabs model with the XSPEC model eye labs (Mihara 
et al. 1990), described by the resonance energy E n s , the res- 
onance width <r res and the resonance depth r res . Using this 
model does not improve the fits significantly, and actually re- 
sults in a slightly worse x 2 /dofs value of 782 / 475 for the 
outburst data. Note that in the case of cyclabs the ratio be- 
tween the energy of the fundamental and first harmonic line 
are fixed to 2, resulting in 1 more degree of freedom. The 
observed E KS of 49. 57$;! keV 1S °f the order of ~ 20% below 
the centroid energy obtained with the gabs model. This dis- 
crepancy stems from the use of a different calculation of the 
line energy and is of the order of 1 0 - 20 % lower than the mea- 
sured gabs energy (see Nakajimaetal. 2010, for details). No 
significant changes in the continuum parameters have been 
observed when using the cyclabs model. A width of cr res = 
19.47f|keV and a depth of r res = 1 -457$;$$ are obtained for 
the outburst observation. In the second observation the res- 
onance energy shows a slight decrease (£ res = 45.77 2 ;okeV), 
and only slightly smaller than the gabs values. These results 
reflect again the difficulty of constraining the possible CRSF 
parameters in the second observation. The width and depth, 
cr res = 5.67$$ keV and r res = 0.467$;}$, are consistent with the 
gabs result showing a decrease in width and depth for both 
observations. 

3.3. Low energy calibration issues 

Strong residuals at lower energies, i.e. between 1 .5 and 2.5 
keV (see Fig. 4) were observed in both data sets. A compar- 
ison with known background properties of the Suzaku/XIS 
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FIG. 4.— Instrumental residuals at lower energies for XIS 1 (blue) and XIS 
03 (red) with the cutoff pi model, a) shows the clearly apparent lines are 
the Au K a line at ~ 2. 1 keV and the Si Kq line at ~ 1.8 keV. Additionally 
a possible Ni K a line at ~ 0.9 which improves the fit slightly, b) shows the 
best fit residuals when the lines are included. 

instrument (Yamaguchi et al. 2006) showed that these fea- 
tures are identical to the known instrumental Si K a and Au 
K q lines, located at 1.74keV and 2.12keV, respectively. In 
many observations of bright sources, e.g. 4U 1907+07 (Rivers 
et al. 2010) and LMCX-3 (Kubota et al. 2010), these energy 
bands are explicitly excluded. In this case, modeling these 
lines with two Gaussian emission features was sufficient to 
minimize the residuals and improve the x 2 /dofs from 1295 / 
483 (without the lines) to the best fit value of 678 / 474. Note 
that the Au K a line at ~ 2.1 keV could be described by the 
same Gaussian emission line at 2.21 ±0.01 keV for both FI 
and BI XIS instruments, whereas the Si K a line appears as 
an emission line for the BI XIS1 at 1.82 ±0.01 keV (Fig. 4, 
blue) and a negative Gaussian line at 1.89 ±0.01 keV for the 
FI XIS 03 combination (Fig. 4, red), indicating that for the 
FI instrument the line is either over-subtracted or not properly 
energy calibrated, leading to a dip in the spectrum. By us- 
ing two independent lines at 1.89keV (FI) and 1.82keV (BI) 
the residuals can be well described. An additional Gaussian 
component at 0.93t$;§2 keV, located close to the Ni K edge at 
0.897keV, slightly improves the residuals (x 2 /dofs of 731 / 
477 without the line, compared to the best-fit value of 678 / 
474) This improvement indicates a possible systematic error 
of the calibration at lower energies for bright sources. This 
line is not very pronounced and can be omitted when using 
the compTT model. 

4. PHASE RESOLVED ANALYSIS 

For a phase resolved analysis the XIS and PIN data of both 
observations were folded with the RXTE determined pulse 
period from Doroshenko et al. (2010) of P sp ; n = 407.719 s, 
Pspin = -4.6 x lO" -7 s/s and the MJD epoch of 54841 .62. Pulse 
profiles with 20 phase bins were created for five different 
energy bands in the 1-10 keV energy range for XIS and 
12-40 ke V energy range for PIN for both observations (Fig . 5 
and Fig. 6). The statistical quality of the GSO data precluded 
the creation of pulse profiles and further spectral analysis. 

During the outburst, the pulse profile consisted of two 
peaks, where the main peak (PI , pulse phase 0.3 -0.6) stayed 
dominant throughout all energy bands and became broader to- 
wards higher energies (Fig. 5). The second peak (P2, pulse 
phase 0.65-0.9) disappeared at energies above ~ lOkeV, 
consistent with the RXTE observations by Doroshenko et al. 
(2010). The third region in Fig. 5 indicates the minimum 
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Fig. 5.— XIS and PIN pulse profiles for different energy bands for the 
outburst observation. PI , P2 and MIN indicate the regions used for the phase 
resolved spectral analysis. 

of the pulse profile (MIN, pulse phase 0.95 - 1 .25), as deter- 
mined from the lower energies of the XIS instrument. 

In comparison, the pulse profile of the second observation 
(Fig. 6) showed a similar Peak PI throughout the whole en- 
ergy band, although narrower at lower energies. No second 
peak was observed at the position of P2, although a small 
'‘bump” in the 5- 12keV energy range was still visible. At 
lower energies, i.e. for the 3 - 5 and 5 - 12keV energy band, 
a small peak showed up on the opposite site of PI , the pulse 
phase where the minimum was defined in the outburst data. 

4.1. Phase resolved spectroscopy 

Spectra were extracted for three different pulse phases 
throughout the outburst: PI , P2 and MIN. XIS and PIN spec- 
tra were used and the same XIS grouping as in the phase av- 
eraged data was applied. The PIN data was again grouped to 
include at least 100 counts per spectral bin. The same spec- 
tral model was used as in the phase averaged analysis and the 
best-fit results are summarized in Table 2. In all pulse phases 
the CRSF component was visible in the spectra. Throughout 
the pulse phase the best-fit values for the galactic and intrinsic 
Ah values did not change significantly. Note that in both mod- 
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FIG. 6.— Same as Fig. 5 for the second observation. The highlighted 
regions indicate the same phase bins as in the outburst observation for a direct 
comparison 

els a slight decrease of the covering fraction during P2 can be 
observed. 

The most significant variation in the cutoff pi model 
was that of the cutoff energy £ cut which decreased from 
29-513 ;|keV in PI, to 1 8.3^5 keV in P2 and down to 
10.71 q 9 keV in MIN. At the same time the power law in- 
dex T varied from 0.42ig;§4 (PI), to 0.56 ^;q? (P2), and to 
— 0. 1 for MIN. In the three phases the CRSF centroid 

energy changed from 66.6^fkeV in PI, and declined to 
55.5^| keV throughout P2, and 52.8!f|keV for MIN. The 
CRSF width declined from 19.21=$ keV in PI to 1 1 .4^ j keV 
in P2, and to 9 . 3 ^j keV in MIN. 

For the compTT model most parameters did not change 
throughout the pulse profile. The observed increase in the 
plasma temperature and the decrease of the optical depth - p is 
a known systematic anti-correlation (see discussion). As with 
the phase averaged data, the CRSF did not show the same 
changes throughout the pulse profile, as in the cutoff pi 
model. The best-fit values for the CRSF centroid eneray were 
57.915;? keV for PI, 57.7!^ keV for P2, and 51.0tfjkeV 


for MIN, and the widths were 13.5ll|^j 5 keV,14.9l3 j keV, and 
9.0l?okeV for PI, P2 and MIN, respectively. Fe K a and 
FeK^ energies were consistent with ~ 6.4 keV and ~ 7.1 keV, 
similar to the cutof fpl values. 

5. DISCUSSION 

This paper presents an analysis of the two Suzaku observa- 
tions of the Be/X-Ray binary 1 A 1 1 18-61 during the peak of 
its outburst in 2009, January and ~ 13 days later. A CRSF, 
detected with RXTE at ~ 55 keV, could be observed in both 
observations, although the significance is lower in the sec- 
ond observation. An Fe K,g line at 7. 13keV has been ob- 
served in addition to the strong, narrow Fe K Q line at 6.4 keV. 
The broad band continuum was modeled with the empirical 
cutof fpl model, including an additional 10 keV systematic 
component to improve the residuals. Softening of the power 
law index F between peak and decay had been observed in 
an earlier outburst and could be confirmed. The Comptoniza- 
tion model compTT, where the lOkeV component was not 
needed, has also been applied. The pulse profiles at lower en- 
ergies changed from a two peaked to a single peaked profile 
between both observations. Phase resolved spectral analysis 
was performed for both observations and the same models as 
in the phase averaged analysis were applied. 

5.1. Outburst behavior 

The third observed outburst of 1A 1118-61 follows a pat- 
tern similar to the second outburst from 1992, i.e. a strong 
peak lasting ~ 3 weeks and an elevated level of emission up 
to 6 weeks afterwards. The time between outbursts was in 
both cases ~ 6200 days, corresponding to ~ 17 years, indi- 
cating that the outburst behavior could be periodic on very 
long time scales. The proposed orbital periods are ~ 350 days 
(Corbet 1986), ~ 58 days (Reig, Fabregat & Coe 1997). and 
most recently 24 days established by Staubert et al. (2010) us- 
ing the delay in pulse arrival time of RXTE monitoring ob- 
servations throughout this outburst. Using the latter method 
the Suzaku light curves are in full agreement with the RXTE 
results (Staubert, priv. comm). A period of 24 days would put 
1A 1 118-61 in the wind accretor region on the “Corbet” dia- 
gram (Corbet 1986), making it a very unique source for a Be 
system. 

A very similar scenario was introduced by Villada et al. 
(1999), who monitored the HJHp emission before and after 
the second outburst and proposed that the optical companion, 
Hen 3-640, has an extended large envelope where a weak 
interaction with the NS can occur. In the scenario of Villada 
et al. (1999), the NS is orbiting the O star in an environment 
with gradually increasing density, until a steady accretion disk 
is created and the X-ray flux suddenly increases. The sudden 
increase in the accretion material would provide the torque on 
the NS to produce the observed changes in the pulse period. 
The surrounding material is then swept out in a short time and 
the system returns to quiescence. According to Villada et al. 
(1999), an interval of 17 years between the outbursts is a rea- 
sonable time scale to accumulate enough material between the 
outbursts. 

5.2. Cyclotron features 

In the outburst, a CRSF has been observed at ~ 55keV 
and the parameters are consistent with the RXTE data 
(Doroshenko et al. 2010). A CRSF in the second observation 
is expected but can not be confirmed with sufficient signifi- 
cance in the established data. 
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TABLE 2 

Phase resolved outburst spectral parameters. Same models as in Table I 


Parameter Cutoffpl Parameter CompTT 

Peak 1 Peak 2 Minimum Peak 1 Peak 2 Minimum 


phabs Ax|[lf) 22 /cm 2 ] 
pcfabs iVH[ 10 22 /cm 2 ] 
covering tract, 
blackbody kT [keV] 
blackbody norm [ 10 ~ 3 ]+ 
£cu. [keV] 

r 

0 66 40 ® 3 
U ‘°°-Q.Q3 

0.52™ 

7 -°3l 

9Q c+7.2 

042^ 

U '^-0.04 

io. 6 !M 

°' 52 AR? 

18.35J 

0.563« 

1 91+0.04 
l - Xi -0.04 

11.734 

0 6 T^ 05 

050^5 

8 T/f-f 

10 7 +1 - 3 
— Q 13-HJ.09 

compTT To [keV] 
compTT kT [keV] 

1 90 +O.O 4 

070^ 

0.543;| 

12.35“ 

1 01 +O.IO 

7.673:1 

1. 22 3 5 

1 1 n+is 

0 59^ 
U-:5V -O.03 

q q+1.6 
1^4 

11. 91/^ 

\ 90+O - 04 
U - DZ -O,02 

10 . Hf 

1 4 a 40 .li 

679^? 

°- /V -0.21 

A. 1 
Pi 

o-issSi 

0-143$ 

o.o 4 5: ( w 

compTT r p 

6 . 313 $ 

4-003$ 

6.603i! 

£crsf [keV] 

66.62SJ 

55.5/p 

1 1 4+0-3 

6 41*^! 

708» 

'• us -0.06 

52.8/1? 

compTT norm 1 

n 1 1 + 0.01 
U - il -0.Q1 
Q+8i7 

13-5^3 4 5 
9 <r fi+4S.4 

0 

7 16+0.04 
'■ 1O -0.05„ 

0-053$ 

57 . 73 ® 

14.9/]® 

7Q e+M.8 
/y,3 -40.2 
6 41+0.01 

70^$ 

^V-O.OS^ 

°.« 

51.0/4® 

<tcrsf [keV] 
fCRSF 

^FeKa [keV] 
£ FcK /3 [keV] 

19 - 2 27 

130.8/™ 

6 41+O.OI 

'• i0 -o.o5 

99 7 +-I 4.6 

7-095 1 , 4 , 

645/85 


Q f)+4.6 
135^.4 

6 42+0-01 

7 19 + 0.14 
/-i “-0.12 
625/85 

Eq. Width K a / [eV] 

373/75 

59/1/11/1 


395/105 

585/13/1 

Flux 2 -iokev[10 36 ergs/sec ]* 

10.5 

8.14 

7.5 





CxiSl/CpiN 

0.98/1.1 

0.97/1.04 

0.99/1.16 


1.04/1.23 

0.99/1.19 

1.08/1.23 

X 2 /dofs 

498/431 

508/420 

474 / 427 


528 / 438 

492 / 423 

542 / 430 


Note. — t Units in Photons keV 1 cm 2 s 1 , $ unabsorbed flux using a distance of 5 kpc 


CRSFs are generated by electrons with energies close to the 
ones determined by the discrete energy states of the Landau 
levels being excited to higher levels followed by de-excitation 
emitting a photon. In this process, photons from the accretion 
column above the magnetic poles with a resonant cyclotron 
line energy, are resonantly scattered out of the line of sight 
creating an absorption line-like feature. These features pro- 
vide a direct measure of the magnetic field strength close to 
the NS surface, where the fundamental cyclotron energy is 
given by 


Lcyc — 


TieB 1 
m e c 1 +z 


1 1 .6 keV 
1+Z 


x -#12, 


(1) 


where Bn is the magnetic field strength near the NS surface in 
units of 10 12 Gauss and z is the gravitational redshift. Assum- 
ing a typical NS mass of 1 .4 M© and NS radius of 10 km gives 
z = 0.3. For 1 A 1 1 18-61 the CRSF at ~ 55 keV corresponds 
to a magnetic field of ~ 6.4 x 10 12 Gauss, making it, together 
with the ~ 50 keV line of A 0535+26, one of the strongest ob- 
served magnetic fields on an accreting neutron star in a binary 
system. 

A possible feature in the residuals at ~ llOkeV could in- 
dicate the existence of a first harmonic line, as observed in 
multiple sources, such as 4U 01 15+63 (Heindl et al. 1999), 
4U 1907+09 (Cusumano et al. 1998; Makishima et al. 1999) 
and VelaX-1 (Makishima et al. 1999; Kreykenbohm et al. 
2002). Including a gabs line at this energy does not signif- 
icantly improve the fit, though a 153 Gd instrumental line at 
~ 150keV, which is due to electron-capture (Kokubun et al. 
2007), as well as possible systematic uncertainties in the back- 
ground weaken the significance of the detection further. 

Assuming that the observed width of the CRSF is due to 
the Doppler broadening of the electrons responsible for the 
resonances, one can estimate the inferred plasma tempera- 
ture in the CRSF resion using equation 4.1 .28 from Meszaros 
(1992): 

/ OkT \ l ' 2 

A+'d = +>Cyc 2 I ! cos ( 0 )l> (3> 

V 2 A+— / 


where Awd corresponds to the CRSF FWHM which is cal- 


culated from <Tcyc. cocyc is equivalent to the CRSF centroid 
energy E Cyc and 0 is the angle between magnetic field lines 
and the line of sight. With cos(0) = 1 the lower limits for the 
plasma temperature are kT « 20.8 keV for the cutoffpl 
model and^T « 12.6keV for the compTT model when using 
the values from Tab. 1 . 

Although the existence of a CRSF in the second observa- 
tion is only marginal, it is very likely that such a line ex- 
ists at lower luminosities. Luminosity dependance of CRSFs 
have been observed in multiple other binary system. In Be- 
Binaries, such as V0332+53, 4U0115+63, and X0331+53, 
where Tsygankov et al. (2006) and Nakajima et al. (2006, 
2010) found a negative correlation between the CRSF cen- 
troid energy and the luminosity of the source for luminosities 
approaching the Eddington luminosity. Staubert et al. (2007) 
on the other hand found a positive CRSF-luminosity corre- 
lation in the low mass X-ray binary HerX-1 for luminosities 
far below the Eddington luminosity. 

The type of correlation seems to depend on whether the ob- 
served luminosity is above or below the Eddington luminos- 
ity. When the luminosity is above the Eddington luminosity, 
the infalling protons start to interact before they are part of 
the accretion column and a “shock front” region, where the 
CRSF most likely occurs, is created. With increasing lumi- 
nosity, the proton interaction occurs farther away from the NS 
surface, where the magnetic field is lower and therefore the 
observed CRSF is seen at lower energies. Below the Edding- 
ton luminosity, the proton interaction does not occur above 
the accretion column but is part of it. When the luminosity 
increases, the accretion “pressure” increases as well and the 
region where the CRSF is created is pressed closer to the NS 
surface, where the magnetic field is higher. This results in 
a positive CRSF-luminosity correlation, as observed in Her 
X-l. 

The flux levels (Tab. 1) obtained for both observations 
showed that the first observation is most likely above, and 
the second observation is definitely below the Eddington lu- 
minosity of the system, so that an anti-correlation between 
CRSF centroid energy and luminosity are expected. 
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5.3. Continuum comparison 

When comparing both observations, changes in the broad- 
band continuum parameters were observed. During the 1974 
outburst Ives, Sanford & Bell Burnell (1975) observed a 
harder spectrum in the peak compared to later observation. 
Looking at the cutoff pi model, the Suzaku data showed 
a similar behavior with a very hard spectrum (power law in- 
dex F = 0.35) at the luminosity peak, and a much softer power 
law index of T ~ 1 in the second observation. Note also that 
the cutoff/folding energy increased with declining luminosity. 
A similar behavior has been observed in a number of differ- 
ent HMXB transients, such as A 0535+26 (Caballero 2009) 
and V 0332+53 (Mowlavi et al. 2006). This is in contrast 
to EXO 2030+375 (Reynolds, Parmar & White 1993), where 
the observed folding energy increased with lower luminos- 
ity. Soong et al. (1990) interpreted this parameter in phase 
resolved results and concluded that the folding energy reflects 
a change of the viewing angle on the accretion column, which 
allows a deeper look into the emission region and therefore 
directly correlates with the observed electron plasma temper- 
ature. The change observed in the phase averaged observa- 
tions follows a similar reasoning, and at lower luminosities 
the emission region is closer to the NS, where the observed 
plasma temperatures are expected to be higher. 

A comparison of compTT parameters between Suzaku and 
RXTE (Doroshenko et al. 2010) showed consistent results 
for the spectral parameters in the outburst. In contrast to the 
RXTE results, the observed column density 1% is higher and 
consistent with the results obtained by the cutof f pi model. 
The main difference between these results and the RXTE data 
lies in the use of a combined Galactic and intrinsic column 
density and the additional need of a low energy black body, 
so that these values cannot be compared directly. Note that 
the RXTE/ PC A instrument used in Doroshenko et al. (2010) 
is not sensitive to data below 3keV and therefore a partial 
covering, as well as a black body component could not be 
modeled. Both, the black body temperature for the low energy 
excess and the photon seed temperature in the compTT model 
decreased by ~ 1/3 between outburst and the second obser- 
vation, indicating a correlation between both parameters. A 
possible explanation is that the source of the compTT seed 
photons is the same as the soft excess modeled by the black 
body component. The observed change in the optical thick- 
ness indicates the the observed material is optically thick in 
the outburst and gets optically thinner for the second observa- 
tion. 

Although an increase in the plasma temperature of the 
compTT model was observed after the outburst, i.e., mainly 
consistent with the interpretation of the cutof fpl param- 
eters, one has to be careful with a direct interpretation of 
the spectral parameter. The optical depth r and the electron 
plasma temperature kT e show a very strong negative correla- 
tion (e.g. Wilms et al. 2006), and the best-fit values cannot 
be used directly for interpretation. Following the definition in 
Rybicki & Lightman (1979), the Compton y parameter: 

4kT c 2 

y= 7 max(r,r ), (3) 

m e c- 

can help in the physical interpretation. Reynolds & Nowak 
(2003) used the Compton v parameter in the description of 
accretion disk coronae around black holes. A value of y « 1 
or slightly higher means that the average emitted photon en- 
ergy increases by an “amplification factor”, A(y) & exp(y), and 
is referred to as “unsaturated inverse Comptonization”. For 


y >• 1 the average photon energy reaches the thermal energy 
of the electrons. This case is called the “saturated inverse 
Comptonization”. Using kT and r p results in y = 2.13 ±0.2 
for the outburst data and y = 1^ 7 for the second observa- 
tion. Typical calculated values of y-parameters are smaller 
than 1, e.g., ~ 0.5 for CygX-1 (Reynolds & Nowak 2003), 
~ 0.2 for 4U 2206+54 (Torrejon et al. 2004), and ~ 0.6 for 
XV 1832-330 (Parmar et al. 2001). Note that the second ob- 
servation is very badly constrained due to the large error bars 
on the electron plasma temperature and the optical depth and 
therefore cannot be used for interpretation. The y parameter 
for the Suzaku outburst indicates that the system is in, or very 
close, to a “saturated inverse Comptonization” state. 

5.4. Fe lines 

In both observations an Fe K a and Fe Tip emission line has 
been observed, where the Fe K Q - Fe K/j normalization ratio 
of 12% is consistent with neutral material. The FeK a EQWs 
were 51 ±4eV and 46 ± 10 eV for the first and second ob- 
servation, respectively, while at the same time the 5-7keV 
power law flux dropped from 6.50 ±0.01 x 10" 10 ergs cm -2 
s" 1 to 1.23 ± 0.0 1± xlO" 10 ergs cm' 2 s -1 . The relatively 
constant EQWs imply that the Fe line emitting region is rel- 
atively close to the source of ionizing flux, so that the ob- 
served Fe line intensity adapts quickly to the changing inci- 
dent flux. The Fe line normalization was 1.79 ±0.06 x 10" 3 
photons cm" 2 s" 1 and 3.15±0.3x 10" 4 photons cm" 2 s' 1 for 
the respective observations. The ratios between Fe normaliza- 
tion and power law flux were consistent in both observations, 
2.75 ±0.10 and 2.56 ±0.27 x 10 6 photons/ergs, reflecting the 
two calculated EQWs. Using the continuum flux difference 
and the time between observations, the calculated change in 
power law flux was ~ 0.415 ±0.001 x 10" 10 ergs cm" 2 s" 1 per 
day in the 12.7 day period between the observations. A linear 
decrease could be assumed using the BAT data of Fig. 1 . For 
the second observation, the Fe normalization was manually 
increased in XSPEC until the resulting EQW matched the up- 
per limit of the measured EQW. This value corresponds to an 
Fe normalization of ~ 3.9 x 10” 4 photons cm" 2 s , assuming 
the measured continuum given above for the second obser- 
vation. Together with the previously calculated constant of 
2.56 x 10 6 , a matching power law flux of ~ 1 .52 could be cal- 
culated. Using this flux difference, together with the rate of 
change in the 5-7 keV energy range, one can then determine a 
maximum time delay which would still preserve the observed 
EQW within errors. The calculated upper limit to the delay 
between the Fe and X-ray emission region is ~ 0.7 days. Note 
that this estimation is very simple and uncertainties from the 
orbital motion or from line of sight assumptions are not taken 
into account here. 

5.5. Phase resolved description 

1 A 1 1 18-61 shows a similar energy and luminosity pulse 
profile dependency as many other sources, e.g. 4U 01 15+63 
(Tsygankov et al. 2007), V 0332+53 (Tsygankov et al. 2006) 
and A 0535+26 (Caballero 2009). In the outburst observation, 
the observed broad double peaked structure at lower energies 
changed to a single peak profile above 10 keV, where the main 
peak (PI) broadened slightly with increased energy and the 
secondary peak (P2) weakened. For the second observation, 
the pulse profile changed significantly. There was still a dom- 
inant primary peak at PI , although it is narrower in the lowest 
energy band. The secondary peak is only marginally indicated 
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TABLE 3 

Pulse fractions for 1 a 1 118-61 


Energy [keV] 

Outburst 

2nd Obs 

1-3 

0.16 

0.73 

3-5 

0,31 

0.90 

5-12 

0.48 

0.53 

12-20 

1.00 

0.78 

20-40 

1.19 

0.83 


between 5-12 keV (see Fig. 6) and absent in the other energy 
bands . On the other hand , between 3-12 ke V a small peak or 
extended shoulder could be observed preceding the primary 
peak. 

To confirm the disappearance of the second peak, pulse 
fractions of the total counts of the P1/P2 regions were cal- 
culated for the indicated energy bands. The pulse fraction 
was defined as (Plcountrate P^countrate) ‘ P2countnttc where the PI 
and P2 regions are indicated in Fig. 5 and 6. Table 3 shows 
that the pulse fraction increases towards higher energies in the 
outburst data, typical for a weakening of the second pulse. In 
the second observation the values do not vary for the differ- 
ent energy bands, and the small peak in the 5- 12keV band 
shows only a marginal smaller pulse fraction than in the other 
energy bands. 

Tsygankov et al. (2007) described similar profiles and pro- 
posed that a misalignment of the rotational and magnetic axes 
of the neutron star leads to the case where one accretion col- 
umn is observed whole, whereas the second accretion column 
is partially screened by the NS surface so that only the softer 
photons, which are created in the higher regions of the ac- 
cretion column, are observed. When the overall luminosity, 
and therefore accretion rate, decreases, the column height de- 
creases and even the soft photons from the second pole are 
shielded. This behavior is in contrast to the variation observed 
in the cyclotron line production region, where the line form- 
ing region is closer to the NS surface at higher luminosities. 

For a more physical picture gravitational effects, such as 
light bending, have to be taken into account when discussing 
pulse profiles (Kraus et al. 1995; Meszaros & Nagel 1985). 
For a canonical NS with a mass of 1 .4 M & and a radius of 
10 km the visible surface of a NS is 83%. With this increased 
surface visibility the pulsed flux from both hot spots is visible 
over a longer part of the pulse phase, and parts of the accre- 
tion column of the second peak are still visible, although that 
hot spot is on the far side of the NS. With decreasing luminos- 
ity, the hot spot size decreases and the visible fraction of the 
second accretion column disappears. 

Pulse profile decomposition methods have been developed 
by Kraus et al. (1995), and have been applied to multi- 
ple sources, e.g. EXO 2030+375 (Sasaki et al. 2010) and 
A0535+26 (Caballero 2011). Under the assumption of a 
slightly distorted magnetic dipole, with this method it is pos- 
sible to disentangle the contribution of the two emission re- 
gions, and constraints on the geometry of the pulsar and on its 
beam pattern can be obtained. The application of this method 
on the 1 A 1 1 1 8 data presented here will be left for future work . 

A more physical approach for a broad band spectrum has 
been introduced by Becker & Wolff (2007), using bulk and 


dynamical Comptonization of photons in the accretion col- 
umn. First results on 4U 01 15+63 have been promising (Fer- 
rigno et al. 2009) and 1A 1118-61 is a good candidate for 
future tests. Schonherr et al. (2007) developed a new physical 
model for the CRSFs based on Monte Carlo simulations for 
the Green’s functions for the radiative transport through a ho- 
mogeneous plasma, which had been successfully applied for 
different sources, e.g., CenX-3 (Suchy et al. 2008). Due to 
advances in the overall code (Schwarm 2007) the test of this 
model on this data set is beyond the scope of this paper. 


6. SUMMARY AND CONCLUSIONS 

In this paper we analyzed the broad band spectrum and 
pulse profiles of 1A 1118-61 during the peak of its third 
observed outburst and compared the results with a second 
observation which occurred ~ 2 weeks after the main peak. 
The time between outbursts is consistent with a continuous 
low level accretion mechanism as suggested by Villada et al. 
(1999). A CRSF during the outburst has been confirmed at 
~ 55 keV indicating one of the highest known B-fields ob- 
served in HMXBs. In the second observation there is only 
a weak indication of a CRSF fundamental line. A change in 
the CRSF centroid energy with respect to luminosity would 
be expected and would help to understand the physical en- 
vironment close to the NS surface. Variations of the CRSF 
can also be observed in the phase resolved spectroscopy of 
the outburst observation. The calculated y parameters show 
that the inverse Comptonization during the outburst is very 
close to saturation. The ratio between the FeK Q and FeKg 
normalization of ~ 12% shows that the emitting material is 
mostly neutral. Using the compTT model, one can deduce 
that the emitting material is optically thick during the out- 
burst and optically thin in the second observation, although 
the known kT-r correlation has to be considered. The pulse 
profile change for different energy bands and luminosities is 
similar to other observed HMXBs and can be explained with 
a misalignment of the rotation axis to the magnetic field axis. 
A change in the pulse profile shape with luminosity has been 
observed and indicates that the visibility of the second hot 
spot changes between the observations. Phase resolved anal- 
ysis throughout the outburst also indicates a change in the 
observed magnetic field, which could be caused by differ- 
ent viewing angles onto the accretion column. Future work 
with a pulse deconvolution technique, as well as a more phys- 
ical model will provide a better understanding of the involved 
physical processes. 
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